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Many eukaryotic cells undergo frequent shape changes (described as amoeboid motion) that enable 
them to move forward. We investigate the effect of confinement on a minimal model of amoeboid 
swimmer. Complex pictures emerge: (i) The swimmer’s nature (i.e., either pusher or puller) can 
be modified by confinement, thus suggesting that this is not an intrinsic property of the swimmer. 

This swimming nature transition stems from intricate internal degrees of freedom of membrane 
deformation, (ii) The swimming speed might increase with increasing confinement before decreasing 
again for stronger confinements, (iii) A straight amoeoboid swimmer’s trajectory in the channel can 
become unstable, and ample lateral excursions of the swimmer prevail. This happens for both 
pusher- and puller-type swimmers. For weak confinement, these excursions are symmetric, while 
they become asymmetric at stronger confinement, whereby the swimmer is located closer to one of 
the two walls. In this study, we combine numerical and theoretical analyses. 

PACS numbers: 47.63.mh, 47.63.Gd, 47.15.G-,47.63.mf 


Some unicellular micro-organisms move on solid sur¬ 
faces or swim in liquids by deforming their body instead 
of using flagella or cilia-this is known as amoeboid mo¬ 
tion. Algae such as Eutreptiella Gymnastica [ij, amoe¬ 
bae such as dictyostelium discoideum^^, but also leu¬ 
cocytes ii and even cancer cells [J use this specific 
way of locomotion. This is a complex movement that 
recently incited several theoretical studies [sj-fij since it 
is intimately linked to cell migration involved in several 
diseases. Some experimental results indicate that adhe¬ 
sion to a solid substratum is not a prerequisite for cells 
such as amoebae 0 to produce an amoeboid movement 
during cell migration and suggest that crawling close to 
a surface and swimming are similar processes. Recently, 
it was shown that integrin (a protein involved in adhe¬ 
sion process) should no longer be viewed as force trans¬ 
ducers during locomotion but as switchable immobilizing 
anchors that slow down cells in the blood stream before 
transmigration. Indeed, leukocytes migrate by swimming 
in the absence of specific adhesive interactions with the 
extracellular environment 131. 


When moving, all micro-organisms are sensitive to 
their environments. Most microswimmers can follow gra¬ 
dients of chemicals (chemotaxis), some microalgae can 
move toward light sources (phototaxis) or orient 
themselves in the gravity field (gravitaxis) [1^, some 
other bacteria move along adhesion gradients (hapto- 
taxis) 01 [l7| . etc. Spatial confinement is another ma¬ 
jor environmental constraint which strongly influences 
the motion of micro-organisms. As a matter of fact, 
in the low-Reynolds number world, amoeboid motion 
generally occurs close to surfaces, in small capillaries 
or in extracellular matrices of biological tissues. Micro¬ 
organisms swim through permeable boundaries, cell walls 
or micro-vasculature. Therefore, the effect of walls on 
motile microorganisms has been a topic of increasingly 
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It has been calculated a long 
and more recently pointed out 


21, 23, 2^, 27, 2^ that swimmers can take advantage of 
walls to increase their motility. Understanding the be¬ 
havior of microswimmers in confinement can also pave 
the way to novel applications in microfluidic devices 
where properly shaped microstructures can interfere with 
swimming bacteria and guide, concentrate, and arrange 
populations of cells [s^. Living microswimmers show a 
large variety of swimming strategies [1^ as do theoretical 
models aiming at describing their dynamics in confine¬ 
ment. 


Felderhof [18| has shown that the speed of Taylor-like 
swimmer increases with confinement. Zhu et al. 2l| used 
the squirmer model to show that (when only tangential 
surface motion is included) the velocity decreases with 
confinement and that a pusher crashes into the wall, a 
puller settles in a straight trajectory, and a neutral swim¬ 
mer navigates. When including normal deformation they 
found an increase of velocity with confinement. Liu et 
al. [ 2 ^ analyzed a helical flagellum in a tube and found 
that except for a small range of tube radii, the swimming 
speed, when the helix rotation rate is fixed, increases 
monotonically as the confinement becomes tighter. Ace- 
moglu et al. [ 2 ^ adopted a similar model but, besides 
the flagellum, they included a head and found a decrease 
of velocity with confinement. Bilbao et aZ.[2^ treated nu¬ 
merically a model inspired by nematode locomotion and 
found that it navigates more efficiently and moves faster 
due to walls. Ledesma et a/. [ 2 ^ reported on a dipolar 
swimmer in a rigid or elastic tube and found a speed 
enhancement due to walls. 


Here, we investigate, by means of numerical and an¬ 
alytical modeling, the effect of confinement on the be¬ 
havior of an amoeboid swimmer, which is a deformable 
object subjected to active forces along its inextensible 
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membrane. Our model swimmer is found to reveal in¬ 
teresting features when confined between two walls, (i) 
We find that straight trajectories might be unstable, in¬ 
dependently of the nature of the swimmer (pusher or 
puller), (ii) For weak confinement, the swimming speed 
can either increase or decrease depending on the confine¬ 
ment strength. For strongly confined regimes, the ve¬ 
locity decreases in all cases recalling previous results on 
different models, (hi) The confined environment is shown 
to induce a transition from one to another type of swim¬ 
mer (i.e. puller or pusher). These behaviors are unique 
to amoeboid swimming (AS) and point to a nontrivial 
dynamics owing to the internal degrees of freedom that 
evolve in response to various constraints. 

The model. Amoeboid swimming is modeled here by 
taking a one-dimensional (ID) closed and inextensible 
membrane, which encloses a two-dimensional (2D) liquid 
of certain viscosity 77 and is suspended in another fluid 
taken to be of the same viscosity, for simplicity. The extra 
computational complexity of dealing with confined geom¬ 
etry restricts our study to 2D, which draws already rich 
behaviors. The effective radius of the swimmer is i?o = 
TT, where Aq is the enclosed area. The swimmer 
has an excess normalized perimeter F = Lo/{2'kRq) — 1 
{Lq is the perimeter) with respect to a circular shape 
(F = 0 corresponds to a circle, whereas large F signifies 
a very deflated, and thus amply deformable, swimmer). 
The strength of conhnement is defined as C = 2Rq/W, 
with W the channel width. Bounding walls are parallel 
to the X direction and y denotes the orthogonal one. 

A set of active forces is distributed on the membrane 
that reacts with tension forces to preserve the local ar- 
clength. The total force density is given by 

dC 

F = Fall - Ccn-b ^t, (1) 

where Fan is the active force to be specified below (which 
we take to point along the normal n for simplicity), C 
is a Lagrange multiplier that enforces local membrane 
incompressibility, c is the curvature, t is the unit tan¬ 
gent vector and s is the arclength. We impose zero 
total force and torque. In its full generality the active 
force can be decomposed into a Fourier series Fa{a,t) = 
Fk{t)e^^°' with a = 27rs/Fo. We first consider 
the case kmax = 3, so that we are left with two complex 
amplitudes F2 and F3. Other configurations of the forces 
have been explored as well (see below). We consider 
cyclic strokes represented by F 2 = F _2 = —Acos{ujt) 
and F 3 = F _3 = Asin(wt), where A is the force ampli¬ 
tude. 

The Stokes equations with boundary conditions (force 
balance condition, continuity of the fluid velocity and 
membrane incompressibility) are solved using either the 
boundary integral method (BIM) or the immersed 
boundary method (IBM) [32|. 

Besides F and C, there is an additional dimensionless 
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FIG. 1: (Color online) Snapshots of an axially moving swim¬ 
mer over time (W = 6Ro). The dashed profiles show a com¬ 
plete period F of deformation and then a few shapes are 
represented over a time of the order of 75Ts. F = 0.085. 


number S = A/itorj), which is the ratio between the time 
scale associated with swimming strokes (Fg = 211 jiS) and 
the time scale of fluid flow due to active force {Tc = y/A). 
Here we take S = 10.0 (the shape has enough time to re¬ 
spond to active forces) and explore the effects of F and C. 
At a large distance from the swimmer, the velocity field 
is governed by Uij = j> FiVjds. Only the (dimensionless) 
stresslet S = — (yyy)/[v/Ts) enters the velocity field 

for symmetric swimmers. E > 0 defines a pusher and 
E < 0 defines a puller. The force distribution defined 
above is found to correspond to a pusher in the absence 
of walls. Below we will see how to monitor a puller or 
pusher and how the walls change the nature of the swim¬ 
mer. 

Results: Axially moving swimmers. We hrst consider 
an axially moving swimmer (AMS) (see Fig. [T]). We 
consider only dimensionless quantities (unless otherwise 
stated). For example, V = VTc/Rq will denote the mag¬ 
nitude of swimming speed. We find an optimal confine¬ 
ment for swimming velocity. Increasing C enhances the 
speed of the swimmer until an optimal (7°, where the 
speed attains a maximum before it decreases. Around 
the optimal value C°, low (high) viscous friction between 
the swimmer and the walls during the forward (recovery) 



FIG. 2: Time-averaged velocity magnitudes (as a function of 
confinement C) of an axially moving swimmer for different F 
values. 
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FIG. 3: Time-averaged (E) as a function of confinement C 
showing the transition from pusher to puller. 



FIG. 5: Time-averaged velocity magnitudes (as a function of 
confinement C) of different swimmers (F = 0.085): migrating 
along one wall (black diamond dashed line), navigating be¬ 
tween two walls (gray circle dotted dashed line) and moving 
along the channel center (black square solid line). The insets 
show characteristic trajectories. 


phase of swimming promotes AMS speed. When the con¬ 
finement is too strong, large amplitude deformations are 
frustrated resulting in a loss of speed. The velocity col¬ 
lapse at strong confinement was also reported for helical 


flagellum [2J, |25| and is expected to occur for all swim¬ 
mer models. Figure [2] shows the swimming velocity mag¬ 
nitude for different F values. That the wall enhances 
motility seems to be a quite general fact, as reported in 


the literature |18l. Il9l. I2ll - l25l| . However, we must stress 


that this is not a systematic tendency. Close inspection 
shows that at weak confinement velocity first decreases 
before increasing, as shown in Ref. 


Results: Swimmer nature evolution. The value of the 
dimensionless stresslet E depends on the instantaneous 
swimmer configuration and its sign instructs us on the na¬ 
ture of the swimmer. We determine the average stresslet 
over a navigation cycle. An interesting result is the ef¬ 
fect of confinement on the pusher/puller nature of the 
swimmer. For small C(< 0.5) the swimmer is found to 
behave as a pusher, while it behaves as a puller for larger 
C(> 0.5). Figure [3] shows the evolution of (E) as a func¬ 
tion of confinement, where a transition from pusher to 



FIG. 4: (Color online) Snapshots of a navigating swimmer 
over time {W = 6Ao). The dashed profiles show a com¬ 
plete period Ta of deformation and then a few shapes are 
represented over a navigation period T of the order of 50Ts. 
F = 0.085. 


puller is observed. 

Results: Instability of the central position. The cen¬ 
tral position after a long time is found to be unsta¬ 
ble. The swimmer exhibits at small C (weakly confined 
regime) a zigzag motion undergoing large amplitude ex¬ 
cursions from one wall to the other. We refer to this as 
a navigating swimmer (NS). Figure [4] shows a snapshot, 
whereas the insets of Fig. [5] display typical trajectories. 
Despite this complex motion, the velocity in Fig. [5] be¬ 
haves with C qualitatively as that of the central swim¬ 
mer. 

The NS trajectory was recently reported 21], in 
the cases of squirmer and three-bead models and also 
observed experimentally for paramecium (ciliated motil- 
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FIG. 6: Average position of the center of mass over a nav¬ 
igation period as a function of confinement C. F = 0.085. 
Circles (diamonds) correspond to the symmetric (asymmet¬ 
ric) motion of the swimmer. The vertical dotted line is the 
demarcation line between a pusher and puller. Note that a 
puller can either navigate or move close to either wall. 













































4 



C 

FIG. 7: Period of navigation as a function of confinement C. 
r = 0.085. Circles (diamonds) correspond to the symmetric 
(asymmetric) motion of the swimmer. 


ity) in a tube 0 pointing to the genericity of naviga¬ 
tion. This instability can be explained analytically (see 
Ref. [Hi). 

A remarkable property is that the navigation mode 
can be adopted both by the pusher and the puller. This 
is in contrast with nonamoeboid motion 2l|, where a 
pusher is found to crash into the wall whereas the puller 
settles into a straight trajectory. These last two be¬ 
haviors are also recovered by our simulations, provided 
the stresslet amplitude is large enough (S^ >> —VDS, 
where D is the dimensionless force quadrupole strength; 
see Ref. (H)- 

Results: Symmetry-breaking bifurcation. At a critical 
C* the symmetric excursion of the swimmer becomes un¬ 
stable and undergoes a bifurcation characterized by the 
loss of the central symmetry in favor of an asymmetric 
excursion in the channel, as shown in the trajectories of 
Fig. [5] (see the supplemental movies in Ref. 33|)- Figure 
shows the average position in y of the center of mass as 
a function of confinement: a bifurcation diagram. This 
bifurcation is very abrupt, albeit it is of supercritical na¬ 
ture. Both slightly before and beyond the bifurcation, 
the swimmer behaves on average as a puller, but still it 
exhibits two very distinct modes of locomotion: naviga¬ 
tion or settling into a quasistraight trajectory (oscillation 
of the center of mass in this regime is fixed by the amoe¬ 
boid cycle). This complexity is triggered by the intricate 
nature of the amoeboid degrees of freedom. 

Results: Other force distributions. Including force 
distributions up to sixth harmonics with various ampli¬ 
tudes leaves the overall picture unchanged, pointing to 
the generic character of AS. The next step has consisted 
in linking the nature of the swimmer to its dynamics. 


We have monitored a pusher or puller type of swimmer. 
If F = 2[sin(a;t) cos(3Qf) — (/3 -|- cos(wt)) cos(2a)]n, we 
have a puller; while if F = 2[(—/3 -|- sin(w<)) cos(3a) — 
cos(wt) cos(2a)]n, we have a pusher (with /3 > 0). (3 
monitors the strength of the swimmer nature (weak and 
strong pusher or puller). We found that for a weak 
enough stresslet, symmetric and asymmetric navigations 
prevail both for pullers and pushers. For a strong enough 
stresslet amplitude (for /3 > /3c ~ 1), we find that the 
pusher crashes into the wall, while a puller settles into a 
straight trajectory. This means that there is a qualita¬ 
tive change of behavior triggered by /3, on which we shall 
report on systematic study in the future. 

Results: Navigation period. The navigation period T 
exhibits a nontrivial behavior with C (Fig. [7]). At small 
C, the period scales as T ~ C~^, and as T ^ at 
intermediate confinement, before attaining a plateau at 
stronger confinement. To dig into the reasons for this 
complex behavior, we provide here some heuristic argu¬ 
ments. In the first regime, the NS swims in a straight 
and monotonous manner towards the next wall. In that 
regime, the period is limited by the distance traveled by 
the swimmer of the order W = 2R/C. This naturally 
yields the C~^ scaling of Fig. [7]for weak C. In the inter¬ 
mediate confinement regime, the magnitude of velocity 
depends linearly on C, so that the period scales as C~‘^ 
(see also Ref. (33|)- After the symmetry-breaking occurs, 
the NS stays close to one of the two walls (inset of Fig. O, 
and its center of mass oscillates with the intrinsic stroke 
period Ts. In this regime, the period is independent of C 
(diamonds in Fig. [7|). 

Analytical results. We have first performed a linear 
stability analysis (H- We find, for small C, that the sta¬ 
bility of the swimmer is governed by the stresslet sign: 
For E > 0 (pusher) the straight trajectory is unsta¬ 
ble, while it is stable otherwise (puller) . For a neutral 
swimmer the trajectory is marginally stable (the stabil¬ 
ity eigenvalue fl, for a perturbation of the form y e^*, 
is purely imaginary). We find that in the intermediate 
C regime the navigation period behaves as C~^. Using 
a systematic multipole expansion the complex behavior 
of the velocity as a function of C (at low C) can be ex¬ 
plained IH - 

Discussion. We believe that the global features re¬ 
vealed by our study will persist in 3D although extending 
our work to 3D simulations will be a challenging task. 
Besides, in order to better match real cells performing 
amoeboid swimming (e.g. leukocytes), cytoskeleton dy¬ 
namics and its relation to force generation will be impor¬ 
tant ingredients to be included in a 3D modeling. 
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